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The pioneering experiments of linear spectroscopy were performed using flames in the 1800s, 
but nonlinear optical measurements had to wait until lasers became available in the twentieth 
century. Because the nonlinear cross section of materials is very small mm, usually macroscopic 
bulk samples and pulsed lasers are used. Numerous efforts have explored coherent nonlinear signal 
generation from individual nanoparticles I3H5] or small atomic ensembles [SHE] with millions of 
atoms. Experiments on a single semiconductor quantum dot have also been reported, albeit with 
a very small yield [9]. Here, we report on coherent nonlinear spectroscopy of a single molecule 
under continuous-wave single-pass illumination, where efficient photon-molecule coupling in a tight 
focus allows switching of a laser beam by less than a handful of pump photons nearly resonant with 
the sharp molecular transition. Aside from their fundamental importance, our results emphasize 
the potential of organic molecules for applications such as quantum information processing, which 
require strong nonlinearities uni. 


Soon after the advent of single molecule spec¬ 
troscopy [HI HU, several groups also reported nonlinear 
measurements on single molecules [TSHIH] . These pio¬ 
neering efforts have detected the incoherent fluorescence 
signal and have not been concerned with optimizing the 
nonlinear yield. However, schemes for many of the future 
quantum technologies rely on optical helds to attenuate 
and amplify other optical signals in a coherent fashion 
and with high efficiency m- In this Letter, we show 
that the intrinsically high efficiency of the linear photon- 
atom coupling also gives access to a highly efficient 
nonlinear interaction of several photons with a quantum 
emitter. 

The schematics of our experiment performed in a he¬ 
lium bath cryostat is depicted in Fig. [T]a and described 
further in the supplementary material. Two laser beams 
at frequencies v-prh and r'pmp were incident on a thin 
naphthalene sample that contained a small concentra¬ 
tion of Dibenzanthanthrene (DBATT) molecules. A solid 
immersion lens (SIL) combined with a simple aspheric 
lens provided a tight focus for the probe beam with 
FWHM^ 320 nm [18]. The probe beam left the SIL with 
a larger divergence so that a part of it was collimated by 
the second aspherical lens. The pump beam, which was 
offset from the lens axis (see Fig. was spatially fil¬ 
tered out by the collecting aspheric lens. In addition, we 
used cross-polarized detection to reduce the contribution 
of the pump to the transmitted probe beam [T5] . 

Figure [^ shows the energy level scheme of DBATT. 
The so-called 0-0 zero-phonon line (OOZPL) with a 
lifetime-limited linewidth of F ~ 20 MHz at the wave¬ 
length of A ~ 619 nm connects the vibrational ground 
states of the electronic ground (S'o) and excited (S'!) 
states [T9|- An important advantage of our system is 
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FIG. 1. Experimental, a. Schematics of the optical arrange¬ 
ment in the cryostat operating at 1.5 K. b. The energy level 
scheme of a DBATT molecule. See text and supplementary 
material for details. 


the large extinction cross section of about 0.5(3A^/27r) 
associated with its OOZPL, where the factor 0.5 signifies 
that about 50% of the population from the excited state 
decays via red-shifted fluorescence channels. 

The interaction of a two-level atom with two opti¬ 
cal fields of amplitudes Ep^p and Aprb at frequencies 
i^pmp and iZpi-b has been described in several works [8l IT^ 
[^ . The bichromatic excitation can be expressed as an 
amplitude-modulated field, resulting in a time-dependent 
Rabi frequency Ut{t) modulated at frequency S = Vp^h — 
^pmp 0. By using a Fourier ansatz for the Bloch vector, 
one can write 
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Poi{t)= ^ (u„ + zr;„)e“‘e2--p-*. 
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where poi denotes the off-diagonal element of the density 
matrix, and Un and Vn are the usual components of the 
Bloch vector. The appearance of the harmonics of 6 in 
Eq. 0 indicates a nonlinear process and exchange of 
energy between the two beams. In an alternative picture 
(see e.g. Fig. and h), the pump field “dresses” the 
molecular levels, giving rise to new transitions at Vp^h = 
z/pnip ± ^^pmp with the pump generalized Rabi frequency 
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FIG. 2. Resonant and near-resonant pump-probe 
spectroscopy, a, b, Transmitted power of a weak probe 
beam as its frequency is scanned across the resonance of a 
single molecule in the absence of a pump beam, c, e, Same 
as in (a) but in the presence of a resonant pump beam at 
two different strengths, d, f, Same as in (b) but in the pres¬ 
ence of a pump beam at A = 5 MHz at two different strengths. 
The zeros on the horizontal axes indicate the pump frequency 
(denoted by the dotted lines). The vertical blue lines show 
the resonance frequency of the unperturbed molecule in each 
case. The pump Rabi frequency is noted in each figure as a 
measure of its intensity. The solid red lines are numerical fits 
using the orders n = —10 to -1-10 in Eq. 0. g, h The dressed 
states of the molecule and the different excitation paths for 
A = 0 (g) and A = 5 MHz (h). The dashed and solid arrows 
depict the probe and pump photons, respectively. The thicker 
and thinner levels in (h) signify the asymmetric distribution 
of populations. Similarly, the darker and lighter arrows rep¬ 
resent more and less probable transitions in each case. 


defined as rip„ip = -I- A^. Here, Hpmp is the Rabi 

frequency caused by the pump and A = Vp^p — Vmo\ 
denotes the detuning between the pump frequency and 
the molecular resonance [2]. 

In Fig. we present transmission spectra of the probe 
laser beam at low Rabi frequency Hpi-b = 5 MHz recorded 
under different pump conditions. The left and right 
columns are for pump detunings of A = 0 and 5 MHz, 
respectively, while Fig. and b show the spectra in the 
absence of the pump beam [18]. The data in Fig. ^ and 
d illustrate that even at a low Rabi frequency well below 
F, there is a clear redistribution of energy from the pump 
to the probe, leading to the development of sidebands. In 
Fig.0 and f these features become clearer, and we ob¬ 
serve a coherent amplification of the probe beam already 
for Hpmp < 2r. 
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FIG. 3. Four-wave mixing, a. The Fourier transform of the 
time-dependent measnrement recorded in transmission. The 
color code nses a logarithmic scale, b, Level scheme of the 
four-wave mixing process. 


In Fig.0, h, we plot the dressed energy levels and the 
different possible pump and probe scattering processes. 
The color coding in Fig. aims to relate the spectral 
features in parts e and f to different pump and probe 
scattering processes in g and h, respectively. The sym¬ 
metric side bands in Fig. 0 are the precursors of the 
Mollow triplet with their zero crossings at in the 

strong pump limit lidHiiin]. In Fig. 0, the feature in 
pink indicates the AC-Stark shift of the extinction dip 
detuned by about Hpj^p, whereas the dispersive feature 
shaded in green represents stimulated Rayleigh scatter¬ 
ing around z/prb ~ t'pmp; a coherent process that involves 
energy transfer between the pump and the probe jT] [53| . 
Finally, at jZprb — ^^pmp + ^^pjj^p a three-photon process 
occurs (highlighted blue), where two pump photons are 
absorbed and one probe photon is emitted in a stimu¬ 
lated process. In this case, the molecule ends up in the 
excited state Si and returns to the ground state after 
a spontaneous decay. Indeed, the resulting fluorescence 
signal was used to detect this hyper-Raman phenomenon 
on a single DBATT molecule for a large pump detun¬ 
ing |16] . Our coherent measurements show that the probe 
is amplified by 0.3%. We note that contrary to the usual 
usage of the dressed picture for the case of strong illu¬ 
mination, Hpmp in our measurements is not much larger 
than F. Consequently, the observed features are less im¬ 
mediately attributable to well-defined transitions among 
dressed states but result from the coherent addition of 
contributions from various channels. 

The interaction of the molecule and the pump photons 
can not only modify the transmitted probe beam, but it 
can also produce light at other frequencies [5] . To investi¬ 
gate this effect, we used a time-correlated single photon 
counting system in a start-stop configuration (see sup¬ 
plementary material) and searched for the beating of the 
probe beam with any other signal. We then calculated 
the Fourier transform of the recorded time-dependent sig¬ 
nal for each pump-probe detuning <5. Figure]^ shows the 
outcome for a measurement at A = 18 MHz. The most 
dominant feature in this figure appearing at frequency 
5 is marked as (i) and stems from the interference of 
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FIG. 4. Time domain measurements, a, Schematics of 
the beating between the pump and probe fields, b-d, Phase 
dependent fluorescence of a single molecule exposed to two 
laser beams with equal Rabi frequencies. The pump beam 
is kept on the molecular resonance while the detuning <5 (see 
legends) of the probe laser is decreased from b) to d). The 
bottom horizontal axes mark time while the upper horizontal 
axes show the phase 5 ■ t of the beating. The blue stripes 
represent the fluorescence lifetime r. The red solid lines are 
fits of our theoretical model to the measured data. The data 
have been corrected for a background linear in the intensity 
of the exciting field caused by off-resonant excitation of other 
molecules. 


the probe beam with residual amount of pump light. In 
addition, however, the data points labeled (ii) reveal the 
existence of a beating signal at 26. Figurej^ sketches the 
four wave mixing process with a degenerate pump beam 
that leads to this new light at frequency r'prb — 2<5 [22] . 

So far we have addressed the steady-state interaction 
between the pump and probe beams as mediated by 
the molecule in frequency space. To visualize the dy¬ 
namic temporal behavior described in Eq. ([^ directly, we 
used the time-dependent measurement described above 
to monitor the population of the excited state by record¬ 
ing the fluorescence signal. To achieve a strong modula¬ 
tion of the excited state population and a large beating 
visibility, we set flprb = f^p^p = 140 MHz. As illustrated 
in Fig. we expect the fluorescence to drop when de¬ 
structive interference reduces the total electric field to 
zero. When the field builds up again, the molecular emis¬ 
sion undergoes Rabi oscillations, which are damped by 
the spontaneous emission to a steady-state level [24] . 
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FIG. 5. Few-photon switching, a. The green data show 
the experimental extinction of an optical beam by a single 
molecule in the absence of a pump field. The blue spec¬ 
trum displays the transmission of the probe beam when the 
molecule was pumped at A = 18 MHz, Hpmp = 50 MHz. The 
red curve projected on the right-hand axis presents the switch¬ 
ing contrast as dehned in the text, b. Same as in (a) but in 
the ideal case of perfect extinction. The green and blue curves 
show the expected spectra for Hpmp = 0, 30 MHz, respectively. 
As displayed by the red curve, the switching contrast reaches 
about 19 dB in this example. The vertical blue lines indicate 
the resonance frequency of the unperturbed molecule while 
the pump frequency is set at zero (denoted by the dotted 
lines). 


Figure j^-d displays three time-resolved fluorescence 
signals, where A was kept at zero and 6 was varied. In 
plots b and c, the Rabi frequency quickly reaches signifi¬ 
cant values within a time much shorter than the fluores¬ 
cence lifetime r = l/27rr, allowing for Rabi oscillations 
to be observed. In both cases, the variation of n{t) with 
time causes a chirp in the temporal oscillations equiv¬ 
alent to subharmonic Rabi resonances in the frequency 
domain [SKIS]. In Fig. [^, the excited-state population 
is not changed for most of the beating period because the 
driving held saturates the molecule within a time com¬ 
parable to r, reducing the visibility of Rabi oscillations. 

Our measurements show that the sharp Fourier-limited 
ZPL of a single molecule can efficiently couple photons 
at two or more frequencies. Such nonlinear interactions 
are of great interest in the context of quantum informa¬ 
tion processing, where single quanta of light and matter 
would be used to transport or mediate information m- 
A key step in many of these endeavors is to switch an 
optical held with only a few, or ideally single, photons. 
We now show that our system readily enters this regime. 
Let us start with the green transmission spectrum in 












































4 


Fig. [^, where a single molecule attenuates a detected 
light signal by 30% (see supplementary material for de¬ 
tails). The blue spectrum in Fig. [§i shows that we can 
recover the transmitted signal at t'moi if we turn on a 
pump field with ripmp = 50 MHz and A = 18 MHz. To 
assess the switching quality, we define Tq^/Toq as the 
switching contrast, where Ton and Tog denote the trans¬ 
mitted probe power with the pump on and off, respec¬ 
tively. As displayed by the red curve in Fig. [^, we 
reach a switching contrast of about 1.8 dB. In the cur¬ 
rent experiment, the molecule required an average of 6- 
30 incident photons within an excited-state lifetime to 
scatter one photon (see supplementary material). How¬ 
ever, theory predicts that this coupling can reach 100% 
in the case of perfect mode matching and unity fluores¬ 
cence branching ratio nzi ng. Figure displays an 
example of high switching performance if we consider a 
probe beam with Hpi-g = F/d, choose A = r/2 and set 


^^pmp = 1.5r. This regime corresponds to a saturation 
parameter + F^/d) ~ 2, implying 

that a pump field with less than a handful of photons 
per T can efficiently switch a beam of light (see supple¬ 
mentary material). 

We have demonstrated the coherent nonlinear interac¬ 
tion of a single organic molecule with only few nearly- 
resonant narrow-band photons, paving the way for single 
photon switching without the need for high-finesse cav¬ 
ities [To]. Our results can also be readily generalized to 
other quantum systems that have a high degree of coher¬ 
ence such as atoms in vacuum [laiiT], rare earth ions in 
crystals |2H| US] or semiconductor quantum dots PEOj. 
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Sample Preparation 

In order to reduce the effect of background fluores¬ 
cence and to optimize the focusing ability of the solid 
immersion lens (SIL), we fabricated our sample in a shal¬ 
low channel etched into a fused silica cover glass (220 nm 
deep and 500 ^m wide) and sealed by the SIL (cubic circo- 
nia). This is done by bringing the SIL and the cover glass 
into optical contact, where they are drawn together by 
van der Waals forces [31] . Subsequently, a small piece of 
naphthalene doped with DBATT was placed on the cover 
glass, next to the SIL and the whole sample was heated 
by a Peltier element until the naphthalene melted. The 
liquid naphthalene was drawn by capillary forces into the 
channel |32j through an opening. The sample crystallized 
after cooling. 


Optical Setup 

The two laser beams were generated from the same 
dye laser with a linewidth of about 2 MHz by using two 
acousto-optical modulators in double-pass configuration 
and variable detuning of —80 MHz to -1-130 MHz. By tun¬ 
ing the laser frequency to the resonances of the molecules 
within the inhomogeneously broadened absorption band 
of DBATT inside a naphthalene crystal, we could address 
each individual molecule separately mm- The trans¬ 
mitted laser intensity as well as the backscattered red- 
shifted fluorescence were measured with single-photon 
counting modules. The pump and probe Rabi frequen¬ 
cies were determined by measuring the broadening of the 
OOZPL in a power-dependent saturation study. 


Spectral Drift Compensation 

The investigated nonlinear effects are very sensitive to 
small frequency variations of only a few MHz. Although 
DBATT in naphthalene is a very stable system, small 
spectral drifts remain. To counter these residual slow 
drifts, the transmission signal was measured during a fast 
sweep of the probe frequency. The spectral position of the 
molecule was measured and the frequency of the pump 
beam was adjusted to maintain the desired relative de¬ 
tuning. This procedure was repeated several 100 times 
until a typical integration time of 1 s per data point was 
reached. 
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Time-resolved measurements 


To register the temporal development of the signal, we 
measured the delay between the arrival of a transmit¬ 
ted photon and a reference signal. The latter was ex¬ 
tracted from the beating of the pump and probe beams 
on an avalanche photodiode in front of the cryostat [55] , 
This was done by using a time-correlated single-photon 
counting system (HydraHarp, PicoQuant) in a start-stop 
configuration. In order to record more than one beating 
period, every fifth stop pulse was taken. The rest of 
the pulses were sorted out by setting an artificial dead 
time using a pulse delay generator (Stanford Research 
DG645). 


Spatially filtered detection 


For the switching experiment, the transmitted light 
was spatially selected by a single mode fiber. Thereby, 
the relative suppression of the pump beam can be im¬ 
proved to more than 200. Additionally, the single mode 
fiber improves the overlap of the probe light and the light 
scattered by the molecule in the detected signal, leading 
to an extinction contrast of about 30%. 


Coupling efficiency 

It is instructive to see what incident photon rate i?inc is 
required to achieve a photon scattering rate i?sca- Com¬ 
bining equations (I) and (3) of Ref. [T7j for steady-state 
monochromatic excitation, one can show that in the case 
of perfect mode matching, 

Rsca, _ ^ 
r? ~ ^ 

-^inc *^0 

7^2 1 

where a = (T+s)‘^o is the scattering cross section 

at a given finite incident power and S = n^/2(A^ -t-T^/d) 
denotes the saturation parameter in the general case of 
detuned excitation [I2l[3l]. Using these results and the 
textbook relation Rsca = Tpee = rS'/2(S' -|- 1) [33], we 
arrive at 

ST/ 4A2\ 

Rinc — -^(^1 + ^j^j- (3) 

Following this derivation, we find that an average of 2 
incident photons per lifetime suffices to reach S = 2 at 
A = F/2 in the ideal case of perfect coupling. 

In our experiment, we determined that 6 — 30 photons 
per excited-state lifetime can suffice as the average rate of 
the incident photons Rinc at the location of the molecule 
for reaching S = 2 on resonance. Here, we carefully cal¬ 
ibrated the optical losses through the cryostat windows 
and lenses on the way to the sample for the case of tight 
focusing. The spread in the measurement stems from the 
variation in the coupling efficiency observed for individ¬ 
ual molecules in the sample. 




